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^—N I The basic concepts of factorizable problems in one-dimensional 

(^ ■ Quantum Mechanics, as well as the theory of Shape Invariant poten- 

. tials are reviewed. The relation of this last theory with a generalization 

Qv I of the classical Factorization Method presented by Infeld and Hull is 

"^ ' analyzed in detail. By the use of some properties of the Riccati equa- 

pi ^1 tion the solutions of Infeld and Hull are greatly generalized in a simple 

■h ; way. PACS numbers: ll.30.Pb, 03.65.Fd. 
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-kj 1 Introduction 

p^ \ The Factorization Method introduced by Schrodinger [iq, 19, 20] and later 



developed by Infeld and Hull |T^ has been shown to be very efficient in 



the search of exactly solvable potentials and his interest has been increasing 
since the introduction by Witten of Supersymmetric Quantum Mechanics 
(SUSY)[^. The bridge between the theory of solvable potentials in one 
dimension and SUSY was estabhshed by Gendenshte'in |^ who introduced 
the concept of a discrete reparametrization invariance, usually called "shape 
invariance" . When studying all these related subjects one is really wondered 
by the almost complete ubiquity of some specific Riccati equations appearing 
in the theory. 



The Riccati equation, which is the simplest first order nonhnear differen- 
tial equation, has a close relation with the group SL{2, M) in the sense estab- 
lished in the celebrated, and unfortunately not so well known as it worths. Lie 
Scheffers theorem [0. This theorem characterizes those first order differen- 
tial equation systems admitting a nonlinear superposition principle. It is also 
known that the problem of finding the general solution for these systems is 
simplified from the knowledge of one or more particular solutions. All these 
aspects have been studied from a group theoretical perspective [|l], Q with 
special emphasis in the Riccati equation [Q, which is nothing but the simplest 
prototype of equation having a nonlinear superposition principle (apart from 
the inhomogeneous linear equations, whose superposition principle reduces 
to a linear one). 

We feel that an appropriate use of the mathematical properties of the 
Riccati equation may be very useful in order to obtain a deep insight in 
the theory of factorizable problems in Quantum Mechanics, as well as in its 
particular class given by Shape Invariant partner Hamiltonians. 

To begin with, the mentioned properties can be used to obtain a simpler 
but more complete presentation, as well as a better understanding of the 
classical results given in |T3. Indeed, we will prove that such results can 



be generalized by simply considering the general solution of certain Riccati 
equation instead of particular ones. In the end, all of the obtained solu- 
tions will give rise to specific, but rather general classes of Shape Invariant 
potentials in the sense of ^. 

Moreover, the techniques to be developed here can be very useful for 
attacking other still unsolved problems. For instance, one could consider 
the study of Shape Invariant potentials depending on several parameters 
transformed by translations as proposed in 0, which is the main subject of 

i- 

The paper is organized as follows: In Section 2 we review the theory of 
related operators and establish the concepts of partner potentials and Shape 
Invariant ones depending on an arbitrary set of parameters. In Section 3 
we establish explicitly the equivalence between a slight generalization of the 
classical Factorization Method [|1^| and the theory of Shape Invariance. Sec- 



tion 4 is devoted to the study of an interesting differential equation system of 
key importance in the development of the subject. The first of its equations 
is a constant coefficients Riccati one for which we will find the solutions in 
full generality. We will use all these results in Section 5, where we obtain 
some rather general classes of factorizable problems which contain as par- 



ticular cases the classical results of |T3. In addition, these will give rise to 



several important families of Shape Invariant potentials which depend on one 
parameter transformed by translation. 

2 Hamiltonians related by first— order diff^er- 
ential operators. The concept of Shape In- 
variance. 

The problem of finding related operators having very similar spectra is now 
a well established subject (see e.g. and references therein). Two linear 
differential operators H and H are said to be A-related if there exists an 
operator A such that AH = HA, where A need not to be invertible. 

Furthermore, if we assume that the operator A relating H and if is a 
first order differential operator, say, of the form 

A=-^ + Wix), (1) 

then, the relation AH = HA, with 



-— + V{x), H = -—^ 
dx^ dx^ 



H = -—,+V{x), H = -— + Vix), (2) 



leads to 

W{V -V) = -W" -V , V -V = -2W' , (3) 

while the relation HA^ = A^H leads to 

w{y ~v) = w" -V' , V -V = -iw . (4) 

One can easily integrate both pair of equations; from the first pair we obtain 
the equation —2WW' = —W" — V and therefore 

V = W'^ -W + c, 

with c being an integrating constant. Following the same pattern with the 
second pair we have 

V = W'^ + W' + d. 



d being also a constant. But taking into account V — V = ~2W' we have 
c = d. We have then the important property that two Hamiltonians H and H 
of the form given by (0) can be related by a first order differential operator A 
given by (|l|) if and only if there exists a real constant d such that W satisfies 
the pair of Riccati equations 

V -d = W^-W, (5) 

V -d = W^ + W. (6) 

Moreover, this means that both Hamiltonians can be factorized as 

H = A^A + d, H = AA^ + d. (7) 

Adding and subtracting equations (^ and (|^) we obtain the equivalent pair 
which relates V and V 

V-d = -{V-d) + 2W'^, (8) 

V = V + 2W'. (9) 

The potentials V and V are usually said to be partners. An important 
concept is the so called Shape Invariance introduced by Gendenshtein p|. He 
supposed that V did depend on a certain set of parameters and considered the 
equations (^ and (^ as a definition of V and V in terms of a superpotential 
W . After, he asked himself what condition was necessary in order to get a 
partner V of the same form as V but for a different choice of the values of the 
parameters involved in V . This relation between V and V is now commonly 
known as Shape Invariance of the potentials 0. 

More exphcitly, we will suppose that our potentials are V = V{x, a) and 
V = V{x,a), where a denotes a set of parameters. Gendenshtein showed 
that if we assume the further relation between V{x, a) and V{x, a) given by 

V{x,a) = V{x,f{a))+R{f{a)), (10) 

where / is an (invertible and different iable) transformation over the set of 
parameters a, then the complete spectrum of the Hamiltonians H and H 
can be found easily. Just writing the a-dependence the equations (|]), (^ 
become 

V{x,a)-d = W^-W, (11) 

V{x,a)-d = W^ + W. (12) 



The simplest way of satisfying these equations is assuming that V{x, a) 
and V{x, a) are obtained from a superpotential function W{x, a) by means 
of 



V{x, a) — d = 


= W^{x,a)-W'{x,a), 


(13) 


V{x, a) — d = 


= W^{x,a) + W'{x,a). 


(14) 



The Shape Invariance property in the sense of [^ requires the further condi- 
tion ([T0| ) to be satisfied. 

Let us remark that the parameter a as well as the transformation law 
/(a) are completely arbitrary up to now, apart from natural requirements as 
differentiability and invertibility. It is clear that the election of a and /(a) is 
what defines the different classes of Shape Invariant potentials. In principle, 
there is no reason why the intersection of these classes should be empty. We 
wiU consider a simple but important type in Section |^. 

3 Equivalence between Shape Invariant po- 
tentials and the Factorization Method 

We consider in this section a slight generalization of the Factorization Method 
as appeared in the celebrated paper |T^. We will prove its equivalence with 



the theory of Shape Invariant partner potentials in the sense of [§ . 

Then, we will deal with the problem of factorizing the linear second-order 
ordinary differential equation 

-^ + r{x,a)y + Xy = 0, (15) 

where the symbol a denotes a set of n independent real parameters, that 
is, a = (ai, . . . , a„). Let us consider a transformation on such parameter 
space /(a) = (/i(a), . . . , /n(a))- We will denote by f'', where fc is a positive 
integer, the composition of / with itself k times. For a negative integer k we 
will consider the composition of f~^ with itself k times and f^ will be the 
identity. The admissible values of the parameters will be f\a), where / is an 
integer restricted to some subset to be precised later. The number A is in 
principle the eigenvalue to be determined. 



In a similar way as in [Q, we will say that (|T3|) can be factorized if it 
can be replaced by each of the two following equations: 

^/-(a)^/-(a)^(^^^) = [X-L{r\a)MX,a), (16) 

HlHly{\,a) = [A - L(a)]i/(A, a) , (17) 

where 

Hl = — + k{x,a), H'^ = -^ + k{x,a). (18) 

ax ax 

Here, k{x, a) is a function to be determined which depends on the set of 
parameters a, and L{a) is a real number for each value of the n-tuple a. The 
fundamental idea of this generalization is expressed in the following 



Theorem 1 Let us suppose that our differential equation jj^ ) can be factor- 



ized in the previously defined sense. Ify{X,a) is one of its solutions then 

yi\,f-\a)) = Hr^''\{\a) (19) 

y(A,/(a)) = Hly{\a) (20) 

are also solutions corresponding to the same A hut to different values of the 
parameter n-tuple a, as it is suggested by the notations. 

Proof.- Multiplying (|16D by hL''^"^ and (P) by H^ we have 

Hf_-Ha) ^rHa) nr^'^^yiX, a) = [\- L{f'\a))]Hr^'^y{\ a) , 

HI HI Hly{\ a) = [X- L{a)]Hly{X, a) . 

Comparison of these equations with ([T6| ) and ([1^) shows that y{X, f~^{a)) 
as defined by (|T^) is a solution of (|T^ with a replaced by f~^{a). Similarly 



y{X, /(a)) given by ( pO]) is a solution with a replaced by /(a). 

It is to be remarked that (K) or (^) may give rise to the zero function; 
actually, we will see that this is necessary at some stage in order to obtain a 
sequence of square-integrable wave functions. 

Indeed we are only interested here in square integrable solutions y{X,a). 
As we are dealing with one-dimensional problems, these solutions can be 
taken as real functions. Under this domain the following Theorem holds: 
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Theorem 2 The linear operators H°^ and H°i are formally mutually adjoint. 
That is, if (pip vanishes at the ends of the interval I , 

f (p{H^ip) dx = f ip{Hl(P) dx . (21) 



Proof.- It is proved directly: 

{H°iip) dx = — I (p —— dx + (p k{x, a) ip dx 
Ji dx Ji 



ip -— dx + (pk(x^a)ip dx = ip(H1(p) dx , 
I dx Jj Jj 

where we have integrated the first term by parts and used that tp(p\si = 0. 

Moreover, it is important to know when (p^Of) and (|20D produce new 
square-integrable functions. 

Theorem 3 Let y{X, a) be a non-vanishing, square-integrable solution of 
(|7^ and ([7^/ The solution y{X, f^^{a)) defined by j^T^ ) is square-integrable 
if and only if X > L{f~^{a)). Similarly, the solution y{X,f{a)) defined by 
l ^m is square-integrable if and only if X > L{a). 

Proof.- It is sufficient to compute 

y{X, f-\a)Y dx = f Hi''^''^y{X,a)Hr'^''^y{X,a)dx 

yiX,a)iHi''^^^HL''^^^yiX,a))dx= {X- Lif-\a))) jy{X,afdx, 



'I J I 

where it has been used Theorem ^ and (|T6D . In a similar way, 

fy{X, f{a)f dx = JHlyiX, a)Hly{X, a) dx 



y{X,a){H^Hly{X,a))dx= {X - L{a)) / y(A,a)'dx, 
I Ji 

where use has been made of Theorem |^ and (0). 

We will consider now the sequence L{f^{a)) and analyze only the cases 
where it is either an increasing or a decreasing sequence. A more complicated 
behavior of L{f^{a)) with k (e.g. oscillatory) will not be treated here. 



Theorem 4 Suppose that L{f^{a)) is a decreasing sequence with no accumu- 
lation points. Then the necessary and sufficient condition for having square- 
integrable solutions of the equations ( |7^ and j ^T^ ) is that there exists a point 



of the parameter space, b = (61, . . . , 6„), such that 

\ = m, HtyiXJib)) = 0, (22) 

provided the function y{L{b), f{b)) so obtained is square-integrable. 

Proof.- Let y{\, a) be a non-vanishing, square-integrable solution of ([T6|) and 
([T7|). In order to avoid a contradiction it is necessary, by Theorem |^, that 
A > L(/~^(a)). If the equality does not hold, one can iterate the process to 
obtain 

y{\ f-\a)f dx={\- L{f-\a))){\ - L{f-\a))) f y{\ af dx . 
I Ji 

Since L{f^{a)) is decreasing with fc, we have that the difference \ — L{f~'^{a)) 
is positive or vanishing and smaller than A— L(/~^(a)). If it still does not van- 
ish, the process can be continued until we arrive at a value /cq such that A = 
L(/-'=»(a)). It is then necessary that |/(A,/-*^° (a)) = i7r*°^"^|/(A, /-'=«+i(a)) = 
0. It suffices to set b = /~'^°(a) to obtain the result. 

Theorem 5 // L[f^[a)) is an increasing sequence with no accumulation 
points, then the necessary and sufficient condition for having square-integrable 
solutions of the equations ( f/^ j and j^PJj) is that there exists a specific point 
of the parameter space, b = {bi, . . . , bn) , such that 

X = m, H'^yiX,b) = 0, (23) 

provided the function y{L{b),b) so obtained is square-integrable. 



Proof.- Let y{X, a) be a non-vanishing, square-integrable solution of ([I6D 
and (0). In order to avoid a contradiction it is necessary by Theorem ^that 
A > L{a). If the equality does not hold, one can iterate the process to obtain 

y(A, f\a)y dx = {X- L(/(a)))(A - L{a)) /y(A, af dx . 
I Ji 



Since L{f^{a)) is an increasing sequence, X — L{f{a)) is positive or vanishing 
and smaller than A — L{a). If it still does not vanish, the process can be con- 
tinued until we arrive at ko such that A = L{f^°~^{a)). Then, it is necessary 
y{\ f°{a)) = Hl''~'^''^y{X, f°-\a)) = 0. It suffices to set b = f°-^{a). 

When L{f^{a)) is a decreasing (resp. increasing) sequence, the functions 
y defined by Hty{L{b), /(&)) = (resp. H'^y{L{b), b) = 0), provided they are 
square-integrable, will be those from where all the others will be constructed. 

We consider now what relation among r{x,a), k{x,a) and L{a) there 
exists. Carrying out explicitly the calculations involved in (0) and (p!?]), 
using (|15D, we find the equations 

k\x, r\a)) + ^^^^'£'^^^^ = -r(x, a) - L{r\a)) , (24) 

fc2(a;,a)-^^^^ = -r(x,a)-L(a). (25) 

We can eliminate r(x, a) between these equations, obtaining 

k\x, r\a)) - k\x, a) + ^^(^'/''(^)) + ^[M^ = L{a) - L{f-\a)) . 

(26) 

Moreover, since (^) and ( pSj) hold for each f^{a), k in the range of integers 
corresponding to square-integrable solutions, we can rewrite them as 

fc^(x,a) + ^^^ = -r{xJ{a))-L{a), (27) 

j^2^^^^^_dk{^ = -r(x,a)-L(a), (28) 

ax 

and from them we can obtain the equivalent pair 

r{x, a) + r(x, /(a)) + 2 k'^{x, a) + 2 L{a) = , (29) 

r{x,a)-r{x,fia))-2^^^^ = 0. (30) 



Both of the equations (p^) and (|25D are necessary conditions to be satisfied 
by k{x,a) and L{a), for a given r{x,a). They are also sufficient since any 
k{x, a) and L{a) satisfying these equations lead unambiguously to a function 
r{x, a) and so to a problem whose factorization is known. It should be noted, 
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however, that there exists the possibihty that equations (^) and ( ^31) did 

not have in general a unique solution for k{x, a) and L{a) for a given r{x, a). 

The equation (EUf) is what one uses in practice in order to obtain results 



of the Factorization Method. We try to solve (|26D instead of ( p4D and (pSf ) 
since is easier to find problems which are factorizable by construction than 
seeing whether certain problem defined by some r{x,a) is factorizable or not. 

Conversely, a solution k{x, a) of (|26| ) gives rise to unique expressions for 
the differences —r{x, /(a)) — L{a) and — r(x, a) — L{a) by means of equations 
([27| ) and (pSD , but it does not determine the quantities r{x,a) and L(a) in 
a unique way. In fact, the method does not determine the function L{a) 
unambiguously but only the difference L{f{a)) — L{a). And this does not 
define L{a) in a unique way at all. To begin with, L{a) is always defined up 
to a constant. And more ambiguity could arise in some cases, as it happens in 
the case studied in [^ . But for the purposes of the application of this method 
to Quantum Mechanics the interesting quantity is L{f{a)) — L{a), as we will 
see below. The same way is undetermined r{x,a), with an ambiguity which 
cancels out exactly with that of L{a) since the differences —r{x, f{a)) — L(a) 
and —r{x, a) — L{a) are completely determined from a given solution k{x, a) 
of®. 

Going back to the problem of finding Shape Invariant potentials in the 
sense of P] which depend on the same set of parameters a, we remember that 



the equations to be satisfied are (|T^) and ( jT^) or the equivalent equations 

V{x,a)-d = -{V{x,a) - d) + 2W^{x,a) , (31) 

V{x,a) = V{x,a) + 2W\x,a), (32) 

as well as the Shape Invariance condition ([T0|). 

Remember that the potentials V{x, a) and V{x, a) define a pair of Hamil- 
tonians 

H{a) = —^ + Vix,a), H{a) = -^ + V{x,a) , (33) 

which can be factorized as 

H{a) = A{a)^A{a) + d, H{a) = A{a)A{a)^ + d , (34) 

where d is a real number and 

d d 

A(a) = — + W{x,a), A\a) = -— + W{x,a) . (35) 

dx dx 

10 



The Shape Invariance condition reads in terms of these Hamiltonians 

H{a) = H{f{a)) + R{f{a)). (36) 

We estabhsh next the identifications between the symbols used in the 
generahzed Factorization Method treated in this section and those used in 
the theory of Shape Invariance. We will see that the equations to be satisfied 
are exactly the same, and that both problems essentially coincide when we 
consider square-integrable solutions. For that purpose is sufficient to identify 

V{x,a)—d = —r{x,f{a)) — L{a), (37) 

V{x,a)—d = —r{x,a) — L{a) , (38) 

W{x,a) = k{x,a), (39) 

R{f{a)) = L{f{a))-L{a), (40) 

and as an immediate consequence, 

A{a) = HI , A\a) = HI , (41) 

for all allowed values of a. In fact, with these identifications it is immediate to 
see that equations (|27|) and (|28|) are equivalent to ([I^ and ([T3|) , respectively. 
Moreover 

V{x,a)-V{xJ{a)) = -r(x, /(a)) - L(a) + r(x, /(a)) + L(/(a)) 

= L(/(a))-L(a)=i?(/(a)), 

which is nothing but equation (p!OD; equations (p9|), (|30D become 

-iV{x, a)-d)- L{a) - (V(x, a) ~ d) - L{a) + 2 W^{x, a) + 2 L{a) 
= -iV{x, a)-d)- {V{x, a)-d) + 2 W'^{x, a) = , 



and 



-{V{x, a)~d)- L{a) + {V{x, a) - d) + L{a) - 2 W'{x, a) 
= -V{x, a) + V{x, a) -2 W'{x, a) = , 



i.e., equations ( pi]) and (|3^), respectively. 

But the identification does not stop here. Let us assume that Theo- 
rem ^ is applicable. We shall see what it means in terms of the Hamilto- 
nians (0). To begin with, we have a certain point of the parameter space 
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b = (61, ...,&„) such that A = L{b) and A^{b)y{L{b), f{b)) = 0, where the 
function y{L{b), f{b)) so defined is square-integrable. We will omit its first 
argument for brevity, writing y{f{b)). It is given by the expression 

yifib)) = iVexp n^ Wi^, b) di\ , (42) 

where A^ is a normalization constant. Note that this wave function has no 
nodes. Since L{f^{a)) is a decreasing sequence, we have that the function 
R{f{b)) = L{f{b)) - L{f-\b)) < for all of the acceptable values of k. 

Then, it is easy to check that y{f{b)) is the ground state of the Hamilto- 
nian H{b), with energy d. In fact, 

H{b)yU\b)) = {A{b)A{by + d)y{f{b)) = dy{f{b)) . 

From equation (|36D we have H{b) = H{f^^{b)) — R{b). The function y{b) is 
the ground state of H{b) with energy d — R{b): 

H{b)y{b) = H{r\b))y{b) - R{b)y{b) = {d - R{b))y{b) . 

Now, the first excited state of H{b) is A{b)y{b): 

H{b)A{b)y{b) = A{b)H{b)y{b) = {d - R{b))A{b)y{b) , 

where it has been used the property H{b)A{b) = A{b)H{b). In a similar way 
it can be proved that A{f~^{b))y{f~~^{b)) is the first excited state of H{b), 
with energy d — R{b) — R{f~^{b)). One can iterate the procedure in order to 
solve completely the eigenvalue problem of the Hamiltonians H{b) and H{b). 
The results are summarized in Table |I|. Note that d has the meaning of the 
reference energy chosen for the Hamiltonians. It is usually taken as zero. 

A similar pattern can be followed when it is applicable the Theorem |^, 
that is, when L{f^{a)) is an increasing sequence. The results are essentially 
the same as when the sequence is decreasing but where now the Hamilto- 
nian with a lower ground state energy is H{b). The basic square-integrable 
eigenfunction y{b) is defined now by A{b)y{b) = 0, that is, 

yib) = Mexp (^- y" Wi^, b) d^^ , (43) 

where M is the normalization constant. Moreover, now R{f^{b)) > for all 
of the acceptable values of k. The results are summarized in Table ^ Again, 
d sets the energy reference level of the Hamiltonians. 
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Eigenfunctions and energies 


H{b) 


H(b) 


Ground state 


yifib)) 
d 


y(b) 

d - R{b) 


A;*'^ excited state 


A{b)---A{f-''+Hb))y{f-''+^(b)) 
d-Er=lRif-''ib)) 


A{rHb))---A(f-^{b))yU'Hb)) 
d-ZLoRif-^'ib)) 



Table 1: Eigenfunctions and eigenvalues of if (6) and H{b) when Theorem^ is 
applicable. The function y{f{h)) is defined by the relation A^{b)y{f{b)) = 0. 



Eigenfunctions and energies 


H{b) 


H{b) 


Ground state 


y{b) 

d 


Vifib)) 

d + R{f{b)) 


A;*'^ excited state 


A^b) ■ ■ ■ AHf-\b))yUHb)) 
d + Y.LiRinb)) 


At(/(6))---At(/fe(b))y{/'=+i(6)) 
d+lZrtlRirib)) 



Table 2: Eigenfunctions and eigenvalues of H{b) and H{b) when is applicable 
Theorem |^. The function y{b) is defined by the relation A{b)y{b) = 0. 
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In both cases the spectra of both Hamiltonians are exactly the same 
(with corresponding eigenfunctions shifted in one step) except for the ground 
state of one of them, which has the lowest possible energy. Only one of 



the eigenfunctions, either (|42|) or (^31) may be square-integrable. It might 



happen, however, that neither of these functions were so. In such a situation 
none of the schemes we have developed would be of use. The conditions on 
the function W{x, h) such that one of the possible ground states exist are 
explained e.g. in 0. Essentially it depends on the asymptotic behavior of 
j"" W{^, h) d^ as X ^ ±00. 

In view of all of these identifications the following result is stated 

Theorem 6 The problem of finding the square integrahle solutions of the 
factorization of (|7^j, given by equations p^j and (jT^j, is the same as to 
solve the discrete eigenvalue problem of the Shape Invariant Hamiltonians 
( p^ in the sense of ^ which depend on the same set of parameters. 

We encourage the reader to compare the results obtained in this section 



with the ones in 0, pp. 24-27], which have inspired this generalization. 

Let us consider now the simplest but particularly important case of having 
only one parameter whose transformation law is a translation, that is, 

/(a) = a — e, or /(a) = a + e, (44) 

where e 7^ 0. In both cases we can normalize the parameter in units of e, 
introducing the new parameter 

m = - , or m = — , (45) 

respectively. In each of these two possibilities the transformation law reads, 
with a slight abuse of the notation /, 

f{m) = m - 1 . (46) 

Then, the equations to be solved for finding Shape Invariant potentials, 
in the sense of [^ , depending on one parameter transformed by a translation 
are 



V{x,m) - 


-d = 


= W^{x,m)-W'{x,m), 


(47) 


V{x,m) - 


-d = 


= W\x,m) + W'{x,m), 


(48) 
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or the equivalent equations 

V{x,m)-d = -{V{x,m) -d) + 2W^{x,m), (49) 

V{x,m) = V{x,m) + 2W'{x,m), (50) 

as well as the Shape Invariance condition 

V{x,m) = V{x,m-l) + R{m-l). (51) 

As a particular case of Theorem ^ we have the following 

Corollary 1 The problem of finding all factorizable problems following the 
Factorization Method stated in / |7^ / is equivalent to find Shape Invariant po- 
tentials in the sense of j^ which depend on one parameter transformed by 
translation. 

The relations among the relevant quantities in both approaches are given 
next for completeness, 



V{x,m) — d = 


= —r(x, m — 1) — L(m) , 


(52) 


V{x,m) — d = 


= —r{x,m) — L{m) , 


(53) 


W{x,m) = 


= k{x,m). 


(54) 


R{m-1) = 


= L{m - I) - L{m) . 


(55) 



We would like to remark that the equivalence between the Factorization 
Method and Shape Invariance has been first pointed out, to our knowledge, 
by several authors almost ten years ago (see e.g. ||2^, ^ and [0). It seems to 
us that most of the authors in these subjects have the feeling (or even a more 
precise knowledge) that such identification exists. But we have not seen so 
far a complete and clear identification in the general case where arbitrary set 
of parameters a and transformation laws /(a) are involved. Our aim is just to 
take a step ahead in the task of clarifying how these methods are interrelated 
since they can be used in more general situations. An important example of 
this is obtained when an arbitrary but finite number of parameters subject 
to translation is involved HI. 
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4 General solution of equations y^ + y' = a, 

zy + z' = h 

In this section we will study the general solution of an ordinary differential 
equation system which will appear as the key point in the solution of the 



problems posed in []T4|, which we will revisit in the next section. 



Let us consider the differential equation system in the variables y and z 

y^ + y' = «, (56) 

yz + z' = b, (57) 

where a and b are real constants and the prime denotes derivative respect 
to X. The equation (|56| ) is a Riccati equation with constant coefficients, 
meanwhile (^TD is an inhomogeneous linear first order differential equation 
for z, provided the function y is known. Recall that the general solution of 
the inhomogeneous linear first order differential equation for v{x) 

-— = a(x)v(x) + bix) , (58) 

ax 

can be obtained by means of the formula 

^,,^,^ rMOexp{-jV,H.He.^ _ (59) 

exp|-j a{^)di\ 



where E is an integration constant. Then, the general solution of (|57|) is 
easily obtained once we know the solutions of (0), i.e. 

expjj y[i)di] 

where we name the integration constant as D. So, let us first pay attention 
to the task of solving (^6]) in its full generahty. 
The general Riccati equation 

— = a2{x)y^ + ai{x)y + aQ{x) , (61) 

where a2(x), ai{x) and aQ{x) are differentiable functions of the independent 
variable x, has very interesting properties. We will recall here some of them 
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which will be of use in our problem. It is a non-linear first order differential 
equation, and in the most general case there is no way of writing the general 
solution by using some quadratures. However, one can integrate it completely 
if some extra information is known. For example, if one particular solution 



yi{x) of (plf) is known, the problem can be reduced to an inhomogeneous 
first order linear equation and the general solution can be found by two 
quadratures. In fact, the change of variable (see e.g. |0, 0) 

u = , with inverse y = yi , (62) 

yi-y u 

transforms (^l]) into the inhomogeneous first order linear equation 

du , . 

— = -[2a2yi + ai)u + a2, (63) 

which can be integrated by two quadratures, for example using (^). An 
alternative change of variable was also found in : 

u = , with mverse y = . (64) 

yi-y u + yi 



This change transforms (|6lD into the inhomogeneous first order linear equa- 
tion 

du [lao \ . . 

-— = \-ai]u + ao, (65) 

dx \ yi J 

which is integrable by two quadratures, as well. We also remark that the 
general Riccati equation ( pTf ) admits the identically vanishing function as a 
solution if and only if ao{x) = for all x. 

Even more interesting is the following property: once three particular 
solutions of (|6TD, yi{x), y2{x), 1/3 (x), are known, the general solution y can 
be written, without making use of any quadrature, by means of the formula 

(y-y^l^y^-y^l = k , (66) 

{y-y2){yi-y3) 

where /c is a constant determining each solution. Solving for y we get 

y2{y3 - yi) k + yi{y2 - ys) .„„v 

y = — 7 w^ • ^^^) 

[ys -yi)k + y2- Vs 
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As an example, it is easy to check that y\k=o = Vi, y\k=i = Us and that the 
solution 7/2 is obtained as the limit of k going to oo. 

The theorem for uniqueness of solutions of differential equations shows 



that the difference between two solutions of the Riccati equation (|6lD has 
a constant sign and therefore the difference between two different solutions 
never vanishes, and the quotients in the previous equations are always well 
defined. 

The equation (^) furnishes a non-linear superposition principle for the 
Riccati equation: there exists a superposition function $(mi,M2,M3, fc) such 
that for any three particular fundamental solutions, the function $(?/i, y2, ?/3, k) 
gives the general solution. 

The first order differential equation systems having this important prop- 
erty are characterized by the so called Lie-Scheffers theorem |T^ , the simplest 



one being the Riccati equation (apart from the inhomogeneous first order lin- 
ear equation, whose superposition principle reduces to a linear one). These 
problems have had a revival after several interesting papers by Winternitz 



and coworkers (see e.g. [^ and references therein), and have been studied 
in |l[] from a group theoretical perspective. In 0] the integrability conditions 
of the Riccati equation, as well as its non-linear superposition principle are 
studied in a unified way by making use of an action on the set of Riccati 
equations. A generalization to other groups and systems admitting such a 
non-linear superposition principle is given in 0. 

We are interested here in the simpler case of the Riccati equation with 
constant coefficients (pBf). The general equation of this type is 

dy 2 , s 

— = a2y + aiy + a^ , (68) 

where 02, ai and oq are now real constants, 02 7^ 0. For a review of some 
of its properties from a geometrical viewpoint see [Q]. This equation, unlike 
the general Riccati equation (^), is always integrable by quadratures, and 
the form of the solutions depends strongly on the sign of the discriminant 
IS. = a\ — 4aoa2. This can be seen by separating the differential equation (|68D 
in the form 

dy dy 

dx . 



a2y^ + aiy + ao , , ai n2 



-^{{y + fkf - i,) 



Integrating (|68|) in this way we obtain non-constant solutions. 
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Looking for constant solutions of (^Hp amounts to solve an algebraic sec- 
ond order equation. So, if A > there will be two different real constant 
solutions. If A = there is only one constant real solution and if A < we 
have no constant real solutions at all. 

We shall illustrate these properties while finding the general solution of 



(56). For this equation the discriminant A is just 4a. Then, the form of the 
solutions depend strongly on the sign of a. If a > we can write a = c^, 
where c > is a real number. The non-constant particular solution 

l/i(x) = ctanh(c(x — A)) , (69) 

where A is an arbitrary integration constant, is readily found by direct inte- 
gration. In addition, there exists two different constant real solutions, 

y2{x) = c, y^(x) = -C. (70) 

Then, we can find out the general solution from these particular solutions 
using the non-linear superposition formula (0), yielding 

B sinh(c(x — A)) — cosh(c(x — A)) . 

B cosh(c(x — A)) — sinh(c(x — A)) 



where B = {2 — k)/k, k being the arbitrary constant in (|67|). Substituting in 
([60| ) we obtain the general solution for z{x), 

^ ^^{B sinh(c(x - A)) - cosh(c(x - A))} + D 
''^'^' 5 cosh(c(x - A)) - sinh(c(x - A)) ' ^ ' 

where D is a new integration constant. 

Let us study now the case with a = in (^6]). By direct integration we 
find the particular solution 

where A is an integration constant. It is clear that now (^) admits the 
identically vanishing solution, and the general solution have to reflect this 
fact. In order to find it is particularly simple the application of the change 
of variable (B3) with yi given by (|73|). Indeed, such a change transforms (|56D 
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with a = into du/dx = 0, which has the general solution u{x) = B, B 
constant. Then, the general solution for (^) with a = is 

y(x) = P , (74) 

with A and B being arbitrary integration constants. If i? = we recover 
the identically vanishing solution as expected. Had we followed the usual 
change of variable (^^ we would have obtained exactly the same result, but 
the calculations would have been a bit longer. Substituting in (^) we obtain 
the general solution for z{x) in this case, 

"^") = l + Bix-A) ' (^') 

where D is a new integration constant. 

The last case to be studied is a < 0. We write then a = — c^, where c > 
is a real number. It is easy to find the non-constant particular solution 

yi{x) = —cta.n{c{x — A)) , (76) 

where A is an arbitrary integration constant, by direct integration. In order 
to find out the general solution, we make the change of variable (|62D or 
alternatively (|64D, with yi{x) given by (^). In both cases the calculations 
are essentially the same and give the general solution of (|56|) for a > 



_ B sin(c(x - A)) + cos(c(x - A)) 
^^^' "" B cos{c{x - A)) - sm{c{x - A)) ' ^ ' 

where B = cF, F an arbitrary constant. Substituting in ( |60D we obtain the 
corresponding general solution for z{x), 

^ \{B sin(c(x - A)) + cos(c(x - A))} + D 
^^"^^ 5 cos(c(x - A)) - sin(c(x - A)) ' ^ ' 

where D is a new integration constant. 

Needless to say, in all of the three cases the solutions can be written in 
many ways, mostly in the cases where exponential, hyperbolic or trigono- 
metric functions are involved. The choice of the form in which the arbitrary 
constants appear might also make the solutions to look a bit different, but 
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these aspects are irrelevant from the mathematical point of view. We have 
tried to give the simplest form for the solutions and in such a way the sym- 
metry between the solutions for the case a > and a < were clearly 
recognized. Indeed, the general solution of (0) for a > can be transformed 
into that of the case a < by means of the formal changes c ^ ic, B ^ iB 
and the identities sinh(ix) = zsin(x), cosh(zx) = cos(a;). The change for B is 
motivated by its definition in the the general solution of (|56|) for a < 0. The 
results are summarized in Table 0. 

We must pay attention to the following point. If we consider, for instance, 
the general solution of ( p6D for a > 0, i.e. equation (|7l|), one could be tempted 
to write it in the form of a logarithmic derivative, 

y{x) = -— log \B cosh(c(x — A)) — sinh(c(a; — A))\ , 
ax 

which is equivalent except for B ^ oo. In fact, if we want to calculate 
lim -— logli? cosh.(c(x — A)) — smh.(c(x — A))\ 

B^oo dx 

we cannot interchange the limit with the derivative, otherwise we would get a 
wrong result. The reason, obviously, is that B cosh(c(x — A)) — sinh(c(a; — A)) 
is not regular as S — > oo. But this limit for B is particularly important since 
when taking it in (fflD , we recover the particular solution ([691) . A similar 
thing happens in the general solutions (|73|) and ([77|), where after taking the 
limit i? — > oo we recover, respectively, the particular solutions ( [73D and (|76D 
from which we have started. Both of ( JT^ ) and ( |77D can be written in the 
form of a logarithmic derivative, but then the limit B —^ oo could not be 
calculated properly. 

The conclusion is the following. If one or more particular solutions of a 
Riccati equation are known, the general solution can be found, for example, 
by one of the methods described above. This general solution depends on 
one parameter characterizing the particular solutions, and in particular one 
should be able to recover the known solutions for some specific values. One 
of these values is usually infinite. If one writes the general solution as a 
logarithmic derivative, the limit when the parameter tends to infinite is to 
be treated with care. 
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Sign of a 


y{x) 


z(x) 


a = c^ > 


B siiih(c(2:-A))-cosh{c(2!-A)) 
'^ B cosh{c(a:-A))-sinh(c(2!-A)) 


|{S sinh(c(a;~A))-cosh(c(a;-A))}+D 
B cosh(c(a;-A))-sinh(c(a:-A)) 


a = 


B 


b{^(x-Af+x-A) + D 
l+B(x~A) 


1+B(x-A) 


a = -c2 < 


, B sin(c(a;-A)) + cos(c(i-A)) 
■^ fl cos{c(i-A))-sm(c(i-A)) 


^{S sin(c(x-A)) + cos{c(a:-A))} + D 
fl cos{c{x- A))-Bin{c{x- A)) 



Table 3: General solutions of the equations ( pq) and (pTl). A, i? and D are 
integration constants. The constant B selects the particular solution of (|56D 
in each case. 

5 The Infeld-Hull Factorization Method re- 
visited: Shape Invariant potentials depend- 
ing on one parameter transformed by trans- 
lation 



We will start this section reviewing the steps of the famous paper [|14| , where 
the Factorization Method was developed in a quite systematic way. It is 
worth mentioning, however, that this method take its roots on previous pa- 
pers by Schrodinger |T^, |T^, ^ and others (see references in [|T^, p. 23]). 
We will apply the mathematical theory developed in the preceding sections 
for solving the problem in a simple way and with full generality, obtaining 
in the end Shape-Invariant potentials in the sense of [Q depending on one 
parameter transformed by translation. 

The key point in the process of finding factorizable problems of type (|T^ 
is to find solutions k{x, a) for the equation (p6D, as we have said in Section |^. 
In our current problem it takes the form 



A;^(x, m + 1) — A;^(x, m) + 



dk{x,m + 1) dk{x,m 



dx 



+ 



dx 



Lirn) — Lim + 1) 



(79) 



which is a differential-difference equation. The idea of solving it in its full 
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generality seems to be very difficult, at least at ffist sight. Instead of doing 
that, it seems to be more sensible to try particular forms of the dependence 
of k{x,m) on x and m. Then, we should find out whether the equation is 
satisfied in each particular case. 

First, note (see [|14[) that there exists a trivial solution of (^), namely 



k{x,m) = f{m) , L{m) = —p{m) , 
where f{m) is any function of m. This gives rise to the problem 



dx"^ 



+ \y = Q. 



which has been discussed completely by Schrodinger |19 



We next try a solution with an affine dependence on m [Q 
k{x,m) = ko{x) + mki{x) , 



(80) 



where ko and ki are functions of x only. Substituting into (^) we obtain the 
equation 



L(m) — Lim + 11 



[(m + l)\kl + k[) + 2(m + l){koki + k'^)] 
-[m^{kl + k[) + 2m(fcofci + k'^)] . 



(81) 



Now we would like to reinterpret the reasoning followed in [|l^ P-27]. Equa- 
tion ( [81| ) reads in its more simplified way 

L{m) - L{m + 1) = 2m{kl + k[) + kj + k[ + 2{koki + k'^) . (82) 

Since L{m) is a function of m alone, the coefficients of the powers of m on the 
right hand side must be constant. Eventually one finds the same coefficients 



to be constant as in the equation appearing after (3.1.4) of |jTj|. Then, the 
equations to be satisfied are 



klko + fcg 



a . 
b. 



(83) 
(84) 



where a and b are in principle real arbitrary constants. When these equations 
are satisfied ( p2D becomes 

L{m) — L(m + 1) = lima + 6) + a . 
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We look for the most general polynomial solution of this equation. It should 
be of degree two in m ii a ^ (degree one if a = 0); otherwise we would find 
that the coefficients of powers greater or equal to three (resp. two) have to 
vanish. Then we put L(r7i) = rm"^ + sm + 1, where r, s, t are constants to be 
determined. Substituting in the previous equation we find the relations 

r = —a, s = —2b, 

and as a result we have the most general polynomial solution for L{m) 

L{m) = —arn^ — 2bm + t , (85) 

where t is an arbitrary real constant. This expression is valid even in the 
case a = 0, being then L{m) = —26m + t. 

In [0, eqs. (3.1.5)] equations (^), ( p4D are written in the slightly more 



restricted way (we use Greek characters to avoid confusion) 



kl + k[ = ~o? 



hko + K = 13, (87) 

where /3 = —70^ if a 7^ 0. This means to consider only negative or zero 
values of a in iW^- Indeed, the solutions of (Bl) for a > are absent in 



T^ , eqs. (3.1.7)], which are supposed to be the most general solutions of the 
system ( |5B| ) and (^7]). However, the solutions appearing when one considers 
the solutions of (^) for a > have their own physical importance. Indeed, 
Infeld and Hull treat particular cases of their general factorization types {A), 



(B) and (E) after having made the formal change a -^ —ia JT^, pp. 27, 30, 
36, 46]. 

But the really important point is that in |T^, even dealing with their 
slightly restricted differential equation system (|^) and (^7]) , they do not 
give the general solutions but simply particular ones, since they only consider 
particular solutions of the Riccati equation with constant coefficients (^61). 
They only consider two such solutions when a ^ and another two when 
a = 0. 

We would like to point out three main aspects now. First, we will treat 



the differential equation system ( ^3]) and (|8^ for all real values of a and 
b. We will find the general solutions of the system by first considering the 
general solution of the Riccati equation (p3|). Second, we will prove that 



all the solutions included in the classic paper |]I4[ are particular cases of 
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that general solutions. Moreover, there is no need of making formal complex 
changes of parameters for obtaining some of the relevant physical solutions, 
since they already appear in the general ones. Thirdly, we will see that rather 
than having four general basic types of factorizable problems [A), [B), (C) 
and {D), where (B), (C) and (D) could be considered as limiting forms of 
(A) [0, p. 28], there exist indeed three general basic types of factorizable 
problems which include the previously mentioned as particular cases, and 
they are classified by the simple distinction of what sign takes a in (0). The 
distinction by the sign of a have indeed a deeper geometrical meaning, but 
we will not go further in this aspect here. See [|, Sec. 4] for more details. 

Moreover, the mentioned lack of generality seems to have been propagated 
to later works trying to generalize the Factorization Method as exposed in 



1^ . See for example some works by Humi |jTO|, [TT|, |T^, [T^. There, more 
general results could be obtained, in principle, by considering negative values 
of certain constant appearing in his reasoning and the general solution of the 
Riccati equation which appears rather than particular ones. For the last two 
of these references, it would be necessary to consider the general solution 
of matrix Riccati equations, which may in turn be formulated by means 
of certain non-linear superposition principle. At this point, it could be of 
practical use part of the extensive work in the field done by Winternitz and 
coworkers (see e.g. [^, ^, |2T| and references therein). 



So, let us find the general solutions of (|8^) and (|J). They are just 
the same as that of the differential equation system ( |5BD and (^7\), simply 
identifying y{x) as ki{x) and z{x) as ko{x), with the same notation for the 
constants. The results are shown in Table §. 

Next we show how these solutions reduce to the ones contained in |T^ . 
For the case a < 0, taking B ^ we recover the factorization type {A) 
of Infeld and Hull |n, eq. (3.1.7a)]. And taking B —>■ i, with a slight 



generalization of the values B can take, we obtain their type (B) (see eq. 
(3.1.76)). For practical cases of physical interest, they use these factorization 
types after making the formal change a -^ —ia [1^, pp. 27, 30, 36, 46]. The 
same results would be obtained if one considers the limiting cases B ^ or 
B ^ 1, respectively, when a > 0, so there is no need of making such formal 
changes. For the case a = 0, taking B ^ oo or B —>■ we recover their 
factorization types (C) and (D) (see their equations (3.1.7c) and (3.1.7(i)), 
respectively. Remember that our convention for the constants appearing in 



equations (^) and (Q) differs slightly from that of equations (3.1.5) of |T^ , 
reproduced here as (|86D and (|87|) with Greek characters for the constants. 
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We show as well some limiting cases of B which give us the particular 
solutions used in the construction of the general ones. Remember that the 
limits B —>■ oo should be taken with care. The arbitrary constant D appear- 
ing in the table is not defined exactly in the same way in all its occurrences 
but it always reflects the fact of having an arbitrary constant wherever it 
appears. 

Let us now try to further generalize (|80|) to higher powers of m. If we try 

k{x, m) = ko(x) + m ki{x) + m^ k2{x) , (88) 

substituting it into ( |79D we obtain 

L{m) - L{m + 1) = Am^kj + 2m'^{3kik2 + 3kl + k'^) 
+2m{kl + 3kik2 + 2kl + 2A;oA;2 + k[ + k'^) + ... , 

where the dots stand for terms not involving m. Since the coefficients of 
powers of m must be constant, from the term in m^ we have k2 = Const. 
From the other terms, if ^2 7^ we obtain that both of ki and /cq have to be 
constant as well. That is, a case of the trivial solution k{x, m) = f{m). The 
same procedure can be used to show that further generalizations to higher 
powers of m give no new solutions [Q. 

Let us try now the simplest generalization of (pOD to inverse powers of m. 
Assuming m 7^ 0, we propose 

k(x,m)= ~^^^^ +ko(x) + mki(x). (89) 

m 

Substituting into ([79|) we obtain 

L{m) - L{m + 1) = - 2 + + . . . , 

m^[m + ly m[m + 1) m[m + 1) 

where the dots denote now the right hand side of (^2]). Then, in addition to 
the equations (|83D and (|8^) the following have to be satisfied 

k\ = e, kok., = f, k'_, = g, (90) 

where the right hand side of these equations are constants. Is easy to prove 
that the only non-trivial new solutions appear when k_i{x) = q, with q non- 
vanishing constant, ko{x) = and ki{x) is not constant. We have to consider 
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Sign of a 


ki{x) and limits 


k()(x) and limits 


Comments 


a = c2 > 


__ B sinh(c(x-A})-cosh(c(x-A)) 
'^ B cosh(c{x - A)) -ainhlclx - A)) 

•^^°°. ctanh{c{a; - ^)) 
— ^^ ccoth{c{x — A)) 


^{S Binb(c(x-A))-coBh(c(x-A))} + D 
B cosh(c(a;~A))-sinh(c(a;-A)) 

B^co^ b ^ -ui 1 A\\ 1 D 


See Q 
See text 
See Q 


"^^°- *> rnthl'rl'r ■11') 1 -^ 


-^=^ ±! + -D cxp{=Fc(x - A)) 


a = 


B 


i,(f (a:-A)2+a;-A) + D 


Type (C) 
Type (D) 


1+B{x-A) 

B^oo 1 

■ x-A 


1 + S(3:-A) 

' 2<.^ '^1 ^ x-A 
-?=^ 6(x A) + D 


a = -c^ < 


B Bin(c{x-A)) + coB(c{x-A)) 
•^ B cos{c{x-A))-sin(c(x~A)) 


^{S sin(c{a:-A)) + cos(c(a:-A))}+D 
B cos(c(a;-A))-sin(c(a;-A)) 


See (JT^) 
Type (A) 
Type (B) 


— ^^ ccot{c{x — A)) 

B^±i .. 
> ±iC 




^^^'. Ti| + Dcxp(Tic{a; - A)) 



Table 4: General solutions of the equations (|83|) and (@), and some limiting 
cases. A and i? are integration constants. The constant B selects the par- 
ticular solution of (^3[) in each case. D is not defined always the same way, 
but always represents an arbitrary constant. 
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Sign of a 


fci (z) and limiting cases 


fco(a;) 


k-i{x) 


Comments 


a = c2 > 


B sinh(c(a:-A))-cosh{c(a:-A)) 
•^ S cosh{c(a:-A))-sinh(c(3:-A)) 






9 eM 


See text 




a = 


S 






9 eK 


Type (F) 


1+B{x-A) 
B^oo 1 


' :c-A 


a = -c^ < 


__ S sin(c(a:-A)) + cos(c(a;-A)) 
■^ S cos(c{x~A))~sinlclx-A)) 

B^O , , , .,. 








Type (£) 





Table 5: New solutions of equations (|83D, ( p4D and (pO|). A is an arbitrary 
constant. B selects the particular solution of (^31) for each sign of a. 



then the general solutions of (^Sj) for each sign of a, shown in Table ^. The 
new results are shown in Table ^. In this table, to obtain really different 
new non-trivial solutions, B should be different from ±1 in the case a > 0, 
and different from in the case a = 0, otherwise we would obtain constant 
particular solutions of (p^). 

For the case a < 0, taking i? ^ we recover the factorization type {E) 
of Infeld and Hull [H, eq. (3.1.7e)]. Again, they use this factorization type 
for particular cases of physical interest after having made the formal change 
a —>■ ia [0, pp. 46, 47]. The same result is achieved by considering the 
limiting case B ^ in a > 0. For the case a = 0, taking B ^ oo we recover 
the factorization type (F) (see their equation (3.1.7/)). For all these solutions 
of (|79D of type (p9D the expression for L{m) is L{m) = —am? — q^ /m? + t, 
with t an arbitrary real constant, which is also valid for the case a = 0. 

It can be checked that further generalizations of ( pOj) to higher negative 
powers of m lead to no new solutions apart from the trivial one and that of 
Tables | and |. 

As a consequence, we have obtained all possible solutions of (^) for 
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k{x,m) if it takes the form of a finite sum of terms involving functions of 
only X times powers oi m. As a consequence of Corollary |l] we have found 
six different, and rather general families of Shape-Invariant potentials in 
the sense of ||^ which depend on only one parameter m transformed by 
translation. These are calculated by means of the formulas (^Tf ), (H), (|5^ 
and (^). We show the final results in Tables ||, and |[ We would like to 
remark here several relations that satisfy the functions defined in Table p. 
In the case a = c^ we have 

/; = c(l - fl) = c{B' - l)hl , h'^ = -C/+/1+ , 

in the case a = 0, 

f, = -Bfl, h', = -Bfoho + l, 

and finally in the case a = — c^, 

/: = c(l + /!) = ciB' + l)h'^ , h'_ = cf^h. , 

where the prime means derivative respect to x and the arguments are the 
same as in the mentioned table, but have been dropped out for simphcity. 

6 Conclusions and outlook 

After a quick review of basic concepts in the theory of factorizable Hamilto- 
nians and Supersymmetric Quantum Mechanics, we have carefully analyzed 
the equivalence between a generalization of the Factorization Method given 
in Ijl^ as to allow the relevant parameters to change in an arbitrary way, and 
the Shape Invariant potentials theory. 

We have treated the particularly simple but important case of only one 
parameter subject to translations, that is, the kind of problems treated by 
Infeld and Hull in their classic paper. To do that, we have considered the 
general solutions of certain Riccati equation with constant coefficients rather 
than particular ones. As a result, we have obtained more general classes of 
factorizable problems (resp. Shape Invariant partner potentials) than the 
ones appearing in ||14|| . 

On the other hand, the bridge beetween Shape Invariance and factoriz- 
able problems has been established more clearly. To this respect, we would 
like to remark that in the interesting paper [^, Sec. VI] a classification of 
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Sign of a 



k(x,m) = ko{x) + mki{x), Lira) 



k{x, m) = q/m + ki{x), L{m) 



a = c^ > 



-f+ix, A, B, c) + Dh+{x, A, B, c) 



-(?rn?- — 2bm + t 



+ mcf+ {x, A, B,c) 



a = 



b ho{x, A, B) + {niB + D)fo(x, A, B) 



'2bm + t 



^+mBfo{x,A,B) 



+ t 



a = -c^ < 



'-f-{x,A,B,c) + Dh^{x,A,B,c) 



crn — 2bm + t 



- mcf—{x, A, B, c) 



c^m^ - -^ + t 



where 

f.f-r A R ^^ — -B sinh{c(x-A))-cosh(c(x-A)) 
j + yx,^,jD,i.j — g cosh(c{x~A))~sinh{c{x-A)) 



f,ix,A,B) = j^^ 



x-A) 



f (rr A n n\ — -B Bin(e(j^-A))+cos(e(3;-A)) 
J^\X,A,a,C) — g coa{c{x-A))-aiTL{c{x-A)) 



h+(x,A,B,c) 
hoix,A,B) = 

h- (X, A, B, C) = g cos{c(a:-A))-siii(c(a:-A)) 



B cosh{c(x — A))—sinh(c(x — A)) 

■{x-A)''+x-A 
1+B(x-A) 



Table 6: General solutions for the two forms of k{x,m) (^OD and (|89D . A, 
B, D, q and t are arbitrary constants. The constant B selects the particular 
solution of (|83D for each sign of a. The constant b is that of (^). 
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Sign of a 


V{x, m) — d, V{x,m) — d, R{m) when k{x,m) = kQ(x) + mki(x) 


a = c2 > 


(t+^af ^2 ^ £^2(6 + ma) + a)f+h+ + (D^ - {B^ - l)(b + ma))hl 
i?(m) = L(m) - L(m + 1) = 2(6 + ma) + a 


a = 


b^hl + (D + mB){D + {m + l)B)/2 + 2b{D + (m + \)B)hho - h 
b^hl + (D + mB){D + {m- l)B)fl + 2b(D + {m ~ \)B)hho + b 
R(m) = L(m) - L(m + 1) = 2b 


a = -c^ < 


_{b+n.af j,2 j^ ^(2(6 + ma) + a)f.h. + {D^ - {B^ + l)(b + ma))hl 
_(b+rp)lf2 j^ ^(2(6 + ^a) - a)f^h^ + (D^ + (B^ + !)(& + ma))hl 
R(m) = L(m) - L(m + 1) = 2(6 + ma) + a 


where f+ = f+{x, A,B,c), fo = Mx,A, B), f- = f.{x, A,B,c) 

h+ = h+{x,A,B,c), ho=ho(x,A,B), h. = h-{x,A,B,c) are defined as in Table | 



Table 7: Shape-Invariant potentials which depend on one parameter m 
transformed by traslation, when k{x,m) is of the form (^0|). A, B, and D 
are arbitrary constants. The constant B selects the particular solution of 
(|83| ) for each sign of a. The constant b is that of (|84D . The Shape Invariance 
condition V{x, m) = V{x, m — 1) + R{m — 1) is satisfied in all cases. 
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Sign of a 



V{x, m) — d, V(x, m) — d, R{m) when k{x, m) = q/m + mki (x) 



a = c^ > 



Xj + m'^c^ + 2qcf+ - m(m + l)c'^{B'^ - l)h\ 
^ + rn^c^ + 2qcf+ - m{m - l)c^{B^ - l)h1 
R{m) = L{m) - L{m + 1) = ^^^ _ ^ + (2m + l)c^ 



a = 



^ + 2gB/o+m(m + l)B2/o^ 
^ + 2qB/o+m(m-l)B2/2 
R{m) = L{m) - L{m + 1) 



(m + l) m'^ 



a = -c2 < 



^ ~ m'^c^ - 2qc/_ + m{m + l)c^{B^ + l)h'i 
^ - m^(? - 2qcf- + m{m - l)c^{B^ + l)h'i 
R(m) = L(m) - L(m + 1) = t-4tt^ - "4- - {2m + l)c'^ 



where /+ = /+(x, A B, c), /o=/o(x,AB), f- = f^{x,A,B,c) 

h^ = h^{x, A, B,c), ho = hQ(x, A, B), h— = h-{x, A,B,c) are defined as in Table I 



Table 8: Shape-Invariant potentials which depend on one parameter m trans- 
formed by traslation, when k{x,m) is of the form (|89|) . A, B, D and q are 



arbitrary constants. The constant B selects the particular solution of ( j83D 
for each sign of a. The constant b is that of (0). The Shape Invariance 
condition V{x, m) = V{x, m — 1) + R{m — 1) is satisfied in all cases. 
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several solutions to the Shape Invariance condition ( pT] ) is given. Comparing 
their ansdstze for the superpotential (6.8) with the one proposed by Infeld 
and Hull, reproduced here as (|89D, is even more clear the relation between 
both approaches. In both of them, the solutions can be generalized simply 
considering the general solutions of a Riccati equation, as we have shown in 
this article. 

But what is even more important is that the use of the properties of 
the Riccati equation provides a great insight in order to study still unsolved 
problems as the one suggested in the end of @, Sec. VI] . That is the subject 
of another article M . 



Finally, we would like to note, since |TJ] is a very referenced and used 
paper, that we have detected one missprint there which may produce later 
unaccurate results. In the expression of the factorization of general Type B 



of [0, page 36], k{x, m) should be d exp{ax) — a{m + c) instead oi d exp{ax) — 
m — c, according to their notation. This missprint is reproduced in their final 
table of factorizations, page 67. However, the function r{x,m) they give for 
that k{x, m) is correct. 
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